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ABSTRACT 


Synthetic Aperture Radar (SAR) data enables direct observa- 
tion of land surface at repetitive intervals and therefore allows 
temporal detection and monitoring of land changes. How- 
ever, the problem of radar automatic change detection 1s made 
more difficult, mainly with the presence of speckle noise. 
This paper presents a new method for SAR image change de- 
tection using the Dezert-Smarandache Theory (DSmT). First, 
a Gamma distribution function is used to characterize glob- 
ally the radar texture data and allows mass assignment throw 
Kullback-Leibler distance. Then, local pixel measurements 
are introduced to refine the mass attribution and take into ac- 
count the context information. Finally, DSmT is carried out 
by comparing the modelling results between temporal images. 
The originality of the proposed method is on the one hand, 
the use of DSmT which achieve a plausible and paradoxi- 
cal reasoning comparing to classical Dempster-Shafer Theory 
(DST). On the other hand, the given approach characterizes 
the radar texture data with a Gamma distribution which al- 
lows a better representation of the speckle. The radar texture 
is being usually modeled by a Gaussian model in previous 
DST and DSmT fusion works. 


Index Terms— SAR image change detection, Kullback- 
Leibler distance, DST and DSmT fusion theories 


1. INTRODUCTION 


Change detection from images covering the same scene and 
taken at different times is of widespread interest for a large 
number of applications, especially in the remote sensing do- 
main. SAR images are very useful tools to surface change 
detection especially in regions where optical data are rarely 
available. In many areas, more and more images are acquired 
on repeated orbits thanks to the repetitivity of radar satellites 
such as ERS 1, ERS 2 and Envisat. However, some difficulties 
are associated with SAR images change analysis: the speckle 
noise, which disturbs automatic change detection. Indeed, the 
speckle imposing a granular texture to radar images makes 
very difficult the use of these images even with slightly differ- 
ent acquisition angles. Nevertheless, despite the complexity 
of data processing, SAR sensors have important properties at 
an operational level, since they are able to acquire data in all 
weather conditions. 


978-1-4244-5654-3/09/$26.00 ©2009 IEEE 


3733 


Recent works have proven that the statistics of SAR im- 
ages can be well modelled by the probability distribution fam- 
ily known as Gamma distribution [1]. Other studies introduce 
some significant statistical measures for the change detection 
purpose [2]. In the same context, an interesting approach 
using statistical change measures between two SAR images 
have been introduced in [3]. The originality of this work is the 
fusion of change signatures computed on two SAR images us- 
ing Dezert-Smarandache theory. This paper extends this work 
and proposes two main contributions. First, SAR texture is 
modelled by a Gamma distribution instead of a Gaussian one. 
Then, in addition to class distribution signatures we introduce 
local change attributes taking into account the context infor- 
mation. 

The paper is organised as follows. Section 2 presents 
some pre-processing tools applied in order to estimate Gamma 
distribution parameters for each class of a temporal image. 
Section 3 exposes global and local measures introduced to 
characterize change signatures. In Section 4, a brief descrip- 
tion of DST and DSmT is presented. The second part of 
this section shows the mass belief assignment used to com- 
pare two images for the change detection purpose. Section 
5 presents experimental results obtained using two temporal 
Envisat images. Conclusions are given in Section 6. 


2. DISTRIBUTION PARAMETERS ESTIMATION 


This section provides an investigation to ascertain the most 
appropriate pre-processing tools which could be used for im- 
age classification in order to estimate Gamma distribution pa- 
rameters. The temporal radar images are first filtered without 
altering SAR texture features and then classified. For each 
class of each image, the Gamma distribution parameters are 
estimated. We can notice here that the filtering process which 
is a delicate pre-processing step is used only for the classifi- 
cation purpose. Parameter estimation is applied on original 
images. 


2.1. Gamma Distribution 


In the case of radar images, previous works showed that the 
Gamma distribution is more accurate than the Gaussian dis- 
tribution to model the real SAR intensity [4]. The Gamma 
distribution PGamma 18 given by: 
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where œ is = shape parameter and ( is the scale parame- 
ter and T(z) = f.~ e™*t7 tdt, z > 0 is the Gamma function. 
Then, a Ax will be characterized by the two Gamma dis- 
tribution parameters qa and /3. 
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2.2. Parameter estimation 
2.2.1. Pre-processing 


Speckle, a form of multiplicative, locally correlated noise, 
plagues imaging applications. For images that contain 
speckle, a goal of enhancement is to remove the speckle 
without destroying important image features. The reducing 
filters have originated mainly from the synthetic aperture 
radar (SAR) community. The most widely cited and applied 
filters in this category include the Lee, Frost, Kuan, and 
Gamma MAP filters. 

More recently, a new Partial Differential Equation (PDE) 
approach to speckle removal was introduced and called 
Speckle Reducing Anisotropic Diffusion (SRAD) [5]. The 
PDE-based speckle removal approach allows the generation 
of an image scale space (a set of filtered images that vary 
from fine to coarse) without bias due to filter window size 
and shape. SRAD not only preserves edges but also enhances 
them by inhibiting diffusion across edges and allowing diffu- 
sion on either side of the edge. 

The aim of SRAD filtering is to ameliorate image classifi- 
cation results in order to perform a better Gamma distribution 
parameters estimation. After SAR image filtering, a K-means 
clustering technique is applied to classify the pixels into k 
and k’ classes for both temporal images. The combination 
between the SRAD filter and the K-means algorithm provides 
a very sufficient classification. 


2.2.2. ML-estimation 


After applying pre-processing steps in order to carry out an 
appropriate classification image mask for both temporal im- 
ages, we estimate statistical Gamma distribution parameters 
for each class. Maximum likelihood estimators are used to 
extract radar texture parameters. 


3. CHANGE MEASURES 


Most of the relevant change detection techniques are based on 
the difference or ratio operators when using radar images. In 
our case, we introduce two measures to characterize change 
signatures. The first one has a global behaviour since it oper- 
ates on classes and the second one is calculated locally for a 
pixel neighbourhood. 


3.1. Gamma Kullback-Leiber distance 


The change detection algorithm is based on the modification 
of the statistics between the two acquisition dates of each 
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pixel. A pixel will be considered as having changed if its 
statistical distribution changes from one image to the other. 
In order to quantify this change, we need a measure which 
maps the estimated statistical distributions for each pixel to a 
scalar index of change. 

In this work, we choose the Gamma Kullback-Leibler di- 
vergence introduced in [6]. The proposed texture similarity 
measure between two Gamma distributions of two image J, 
and Tə can then be given by: 


Dict Gonna li Iz) = log D(a) T 4) log BU») 
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For simplicity, the distribution value J; (i, 7) for the pixel 
(i, j) of the image Iņ, k = 1, 2, is noted J;, in the above equa- 
tion. The factors a/*) and 3), are the shape and the scale 
parameters associated with distribution J;,, k = 1, 2. 
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3.2. Local pixel measures 


The contrast measure takes care of the pixel realizations from 
a stochastic point of view. To highlight contrast information, 
the conventional Pixel by Pixel Ratio Measure (PPRM) be- 
tween two SAR images is exploited. This detector is well- 
known and widely used in SAR imagery due to its ability to 
greatly reduce the speckle influence on the change map. The 
PPRM of a pixel between two image J; and [> is computed 
and normalized on a small window through a contrast mea- 
sure r(i, j) defined by: 
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r(i, j) = log(max( a a: 


where v is the pixel (i, j) neighborhood defined by a given 
window. 

Moreover, we define a significance measure as a criteria 
dedicated to reduce the ambiguity of the pixel behaviour. The 
significance measure of a pixel (i, j) for images J, and Ib 
may be defined by: 
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We can notice that the contrast and the significance measures 
are calculated for a pixel neighbourhood and not for a class 
distribution. Besides, those measures characterize the evo- 
lution of each pixel and are not necessary linked to ground 
change but give potential candidates to ground evolution. 
That is why those indicators are to be considered into an 
evidential and paradoxical reasoning. 


4. THE DEZERT-SMARANDACHE THEORY 


The Dezert-Smarandache Theory (DSmT) [7] is a general- 
ization of the classical Dempster-Shafer Theory (DST) [8], 
which allows formal combining of rational, uncertain and 
paradoxical sources. The DSmT is able to solve complex fu- 
sion problems where the DST usually fails, especially when 
conflicts between sources become large. In this section, we 
will first review the principle of the DST before discussing 
the fundamental aspects of the DSmT. 


4.1. Dempster-Shafer Theory principle 


The DST makes inferences from incomplete and uncertain 
knowledge by combining sources of confidence, even in 
the process of partially contradictory sensors. In the DST, 
there is a fixed set of mutually exclusive and exhaustive 
elements, called the frame of discernment, which is sym- 
bolized by the set of N elements potentially overlapped 
© = {61,62,...,9n}. The frame of discernment © defines 
the propositions noted A; for which the sources can provide 
confidence. Information sources can distribute mass values 
on subsets of the frame of discernment, A; € 2°. If an 
information source can not distinguish between two proposi- 
tions A; and A,, it assigns a mass value to the set including 
both hypotheses (A; U A;). The mass distribution for all 
hypotheses has to fulfill the following conditions 


it. ma) =0. 


iii. Toy c20 M(Ai) = 1. 


Mass distributions from d different sources are combined with 
Dempster’s orthogonal rule. 
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where K = dian. nay ME, m;( Ai) is a normaliza- 
tion factor measuring the conflict between the sources. 
Two functions can be evaluated to characterize the uncer- 
tainty about the hypotheses. The belief function Bel( A) = 
> a,ca M(Ai) measures the minimum uncertainty about A, 
whereas, the plausibility function PI(A) = > 7 4,,449 M(Ai) 
reflects the maximum uncertainty value. 


4.2. The Dezert-Smarandache Theory (DSmT) 


While the DST considers © as a set of exclusive elements, the 
DSmT relaxes this condition and allows for overlapping and 
intersecting hypotheses. This allows for quantifying the con- 
flict that might arise between the different sources throughout 
the assignment of non-null confidence values to the intersec- 
tion of distinct hypotheses. 


4.2.1. Definition 


Let © = {61,02,...,9n} be a set of N elements which 
can potentially overlap. We consider in our example © = 
{Ocn;9no.ch} involving only two elementary hypotheses 
*change’and ’no change’. The hyper-power set DÈ is de- 
fined as the set of all composite hypotheses obtained from © 
with N and U operators such that: 


is 0, 01,00, ..., 0N eD’ 
ii. If A,B € D? then (AN B) € D? and (AU B) € D®. 


iii. No other elements belong to DÌ except those defined in 
2. and ti. 


We define D® in our example by DÌ = Vs Ochs Oneccns een) 
Ome chi Och U Oroshi: 

As in the DST, the DSmT defines a map m(.) : D° — [0,1]. 
This map defines the confidence level that each sensor asso- 
ciates with the element of DÌ. This map supports paradoxical 
information and we have this condition: 


XO m(A) =1. (6) 
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The belief and plausibility functions are defined in the same 
way as for the DST. The DSmT rule of combination of con- 
flicting and uncertain sources is given by the above equation: 


d 
Ss  [ mi(4i). (7) 
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4.2.2. Mass belief assignment 


Changes may be expected for a pixel (7, 7) when the Gamma 
Kullback-Leiber distance Dg tGamma(, I2), the contrast 
measure r(i, j) and the significance measure s(i, j) become 
significant. Furthermore, the decision may be taken when 
there is a contradiction between the three indicators. Then, 
we propose in Table 1 the mass assignment which is done 
by an appropriate combination of change signatures. The 
decision is taken through the maximum of credibility with 
Belief on change Bel (Oen). 


Hypothesis Mass belief assignment 
r(i, j) x DriGamma dis Iz) xX s(i, j) 


(1 — (i,j) * (1 — DkLGamma (Iı, I2)) * (1 = s(i, 9) 
r(i, j) * (1 — DrkLGamma(lı, T2)) * s(i, j) 
1— (m (Ø) +m (Och) +m (Ono-ch) +m (Och U Onoch)) 


Table 1. Mass Belief Assignement for image fusion. 
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Fig. 1. Original Envisat SAR images. 


5. RESULTS 


The proposed method has been applied on two temporal 
Envisat SAR images covering the same region of Tunis City, 
North Africa (see Fig. 1). The first one is acquired in 2005 
whereas the second one is acquired in 2007. A considerable 
amount of urban changes occurred between these two dates. 
These two images are already registered. 

Fig. 2 shows a comparison between the change measures 
results obtained respectively with the KL distance, the con- 
trast and the significance measures. As we can see, some 
change regions are highlighted by the three signatures. How- 
ever, others are enhanced only by one or two measures which 
underline the complementarities of these attributes. 

Fig. 3 shows decision results obtained by maximizing 
the belief change function exposed previously. Even if this 
method takes a strong decision, it 1s possible to analyse the 
belief response or the interval [credibility, plausibility] to in- 
troduce a confident interval into the decision. This decision is 
made under DST and DSmT. 

Changes detected by DSmT processing (see Fig. 3 (Mid- 
dle) shows that this approach provides better performances 
than the classical DST one (see Fig. 3 (Left)). Comparing 
these results to original images presented in Fig. 1, we can 
notice that the proposed method detects the most significant 
changes with a very few false alarms. So, taking into account 
conflict between the different sources enhances considerably 
the change detection results. Fig. 3 (Right) gives a quanti- 
tative evaluation (ROC plots) of the results using the ground 
truth. It shows the performances of all used measures and 
it confirms that DSmT (blue) is more accurate comparing to 
DST (red) and all the other indicators. 


6. CONCLUSION 


Change detection in multi-temporal SAR images requires the 
computation of specific measures which have to be sensi- 
tive to various kinds of changes and robust to speckle effects. 
Those measures are combined and used by DSmT to take into 
account conflicts between sources which enhance change de- 
tection results. The performance of the proposed technique 
is very suffcient comparing to DST results and original im- 
ages. In future works, we can take the benefit from the DSmT 
hybrid model and introduce the multiscale information in the 
proposed change detection approach. This will allow detec- 
tion of changes with various sizes. 


Fig. 2. Change measures results given respectively (from left 
to right) by Kullback-Leiber distance, contrast and signifi- 
cance measure. 

















Fig. 3. Change detection results. (Left) DST Changes (Mid- 
dle) DSmT Changes (Right) ROC plots comparison between 
Contrast, signifiance, KLGamma, DST and DSmT change de- 
tection measures. 
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